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Some of the most intriguing problems in solid state physics arise when the motion of one electron 
dramatically affects the motion of surrounding electrons. Traditionally, such highly-correlated elec- 
tron systems have been studied mainly in materials with complex transition metal chemistry. Over 
the past decade, researchers have learned to confine one or a few electrons within a nanoscale semi- 
conductor "artificial atom", and to understand and control this simple system in exquisite detail. 
• In this Article, we combine such individually well-understood components to create a novel highly- 

£\j ' correlated electron system within a nano-engineered semiconductor structure. We tune the system 

in situ through a quantum phase transition between two distinct states, one familiar and one subtly 
new. The boundary between these states is a quantum critical point: the exotic and previously 
elusive two-channel Kondo state, in which electrons in two reservoirs are entangled through their 
interaction with a single localized spin. 

The Kondo effect has become a hallmark of coherent electron transport in a variety of nanostructures ranging 
from lithographically-defined semiconductors Q to carbon nanotubes and molecules 0, . Kondo applied the 
phenomenological Hamiltonian Q 
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Hk — Js - S + H reservoir (1) 



' to describe a magnetic impurity embedded in a host sea of electrons. A localized spin S couples antiferromagnetically 
c*j ■ with strength J to spins s of electrons in the surrounding reservoir. H reS ervoir represents the free electrons in the 
reservoir. At temperatures below the Kondo temperature Tkj electrons in the reservoir screen the localized spin. The 
Kondo Hamiltonian was later found to be derivable from the more microscopic Anderson model, which consists of 
an electron bound to an impurity site in a metal host (Fig.^a)). Here, Kondo's antiferromagnetic coupling emerges 
t-H , from tunneling on and off the local site. 

Many systems of strongly-interacting particles can be understood in the framework of Landau's Fermi liquid the- 
ory, whose basic entities, termed quasiparticles, roll most effects of interactions into changes in particle properties 

^ i 1 such as mass and energy. Although the Kondo ground state is complex, its excitations can still be described as 

weakly-interacting quasiparticles. Some of the most intriguing problems in solid state physics arise when this simpli- 
t-H ' fication cannot be applied. Examples of such highly-correlated systems include Luttinger liquids, fractional quantum 
Hall Laughlin liquids, high-temperature superconductors, and the two-channel Kondo system, a novel state studied 
experimentally in this Article. 

In the two-channel Kondo (2CK) model, introduced 25 years ago by Nozieres and Blandin, and independently 
by Zawadowski 0, , a localized spin S is antiferromagnetically coupled to two independent reservoirs of electrons 
I ■ according to the Hamiltonian 

-o 

£h [ H2CK — J\S\ ■ S + J2S2 ' S + ^reservoirs- (2) 

q , The symmetric 2CK state is formed when the two independent channels (or reservoirs) are equally coupled to the 
magnetic impurity, i.e. J\ = J 2 . Each reservoir individually attempts to screen the local spin. Since they cannot 
both screen the spin, this is an unstable situation, resulting in a new ground state in which the local impurity is 
only partially screened. Unlike for single-channel Kondo (1CK), in the 2CK state the quasi-particle concept of Fermi 
liquid theory does not apply: the decay rate for a low energy excitation (~ y/e) is greater than the energy e of the 
excitation itself. Stable low- lying excitations must thus be collective 0- 

Any difference in channel coupling will force the system away from the non-Fermi liquid 2CK state and toward the 
1CK state associated with the more strongly-coupled reservoir. Although the symmetric 2CK state might therefore 
seem difficult to access, it has been invoked to explain remarkable low-energy properties of some heavy Fermion 
materials^ [H E3 and glassy metals[H Q El However, the connections of these experimental observations to 
2CK theory remain problematical , in part because the microscopic electronic structure of the various materials is 
unclear. 

In this Article, we present experimental results on an artificial impurity which is designed to display 2CK. Crucially, 
we can precisely control the coupling constants J\ and J2, while maintaining the independence of the two channels. 
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FIG. 1: Single Channel Kondo Effect (1CK). (a) The Anderson model describes a magnetic impurity in a metal as a single spin-degenerate 
state (right side of green barrier) coupled to a Fermi reservoir of electrons (left). Coulomb interaction U between localized electrons favors 
having only a single electron in the localized state. The antiferromagnetic coupling J betwe en the localized spin and the reservoir depends 
on the tunneling rate T, the depth of the level eo, and U, according to J ~ TU / (eo( e O + E/D I21I . At low temperature, high-order tunneling 
processes between the local state and the Fermi reservoirs coherently add together to screen the localized electron spin, (b) Two Channel 
Kondo Effect (2CK). A localized electron is now coupled to two independent Fermi reservoirs (blue and red). If the two independent 
reservoirs are equally coupled to the localized spin, each will individually attempt to screen the spin, resulting in the formation of a 
highly-correlated electron state, (c) Physically separating two reservoirs does not suffice to make them independent. If a localized electron 
can hop off the site to the right reservoir and a new electron can hop onto the site from the left, the two reservoirs will cooperate in 
screening the localized spin. To create two independent screening channels, processes which transfer electrons from one reservoir to another 
must be suppressed, (d) Experimental Realization of 2CK. We add an additional finite reservoir (red) to an artificial magnetic impurity 
connected to an infinite reservoir comprised of two conventional leads (blue), (e) Coulomb blockade suppresses exchange of electrons 
between the finite reservoir and the normal leads at low temperature. The two reservoirs (blue and red) hence act as two independent 
screening channels (see Text.) 



The system is built from a GaAs/AlGaAs heterostructure containing a low density (n e = 2 x 10 11 e~ /cm 2 ), high 
mobility (/x = 2 x 10 6 cm 2 /Vs) two-dimensional electron gas (2DEG) 68 nm below the surface. We follow the proposal 
by Oreg and Goldhaber-Gordon^^ to produce two independent screening channels for an artificial magnetic impurity. 
A gate-defined quantum dot containing ~ 25 electrons in an area of 0.04 /im 2 plays the role of our magnetic impurity 
(Fig. ^d), left). Its bare charging energy U w 1 meV and its average single-particle level spacing A w 100 \xeV . 
Previous experiments which claim to probe 2CK pH[l^ll^ll^ll4lll5T | used a local orbital degeneracy in place of spin, 
freeing s pin of the surrounding conduction electrons to act as the channel index. In contrast, our local degeneracy is a 
real spin[l3, and we use two physically-separated reservoirs (red and blue in Fig.[IJd)) for the screening channels. Two 
leads of the small quantum dot cooperate as a single screening channel with antiferromagnetic coupling J lT ("infinite 
reservoir" ) [lflj . An additional lead is made finite in size (red, Fig.QJd)), so that adding or removing an electron from 
this reservoir is energetically forbidden at low temperature, a phenomenon known as Coulomb blockade. The area of 
the finite reservoir is ~ corresponding to a charging energy E c = 100 [ieV w 1.2 K, and a single-particle level 

spacing Af r = 2 fieV w 25 mK ("finite reservoir"). This level spacing is only slightly larger than the base electron 
temperature of our dilution refrigerator - 12 mK, as determined by Coulomb blockade thermometry on the small 
quantum dot - and indeed we cannot resolve these levels even at base temperature. Hence, the finite reservoir has an 
effectively continuous density of states and can screen the magnetic impurity [2(lj . Since Coulomb blockade prevents 
exchange of electrons with the other leads, the finite reservoir acts as a second Kondo screening channel (Fig. [He)), 
with antiferromagnetic coupling Jf r , allowing the possibility of observing and studying 2CK. The 2CK Hamiltonian 
(Eq. (0)) has three possible ground states, depending on the relative couplings to the two reservoirs: 1CK with the 
finite reservoir ( Jf r > J; r ), 1CK with the infinite reservoir (Jj r > Jf r ), and 2CK at the quantum critical point Jf r = J lr . 

In Fig. we demonstrate that the small quantum dot can act as a tunable magnetic impurity and display the 
single-channel Kondo effect. If the small quantum dot has an odd number of electrons, it has a net spin and acts as a 
magnetic impurity. With gate n de-energized (0 V), the system has three conventional leads (blue and red in Fig-QJ, 
all of which cooperate to screen the magnetic impurity with a single energy scale folk- At temperature T < Tk the 
Kondo effect enhances scattering and hence conductance from one lead to another. We measure the conductance 
g = d//(iVds|vd B =o between the two blue leads. As temperature is increased, the Kondo state is partially destroyed, so 
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FIG. 2: (a) Scanning electron micrograph of a device similar to that measured. The device consists of a small quantum dot (magnetic 
impurity, left) coupled to conventional leads (top and bottom left) and to a large quantum dot (finite reservoir, right). Electrons are 
depleted under every gate by application of negative voltages. In the experiments described here, voltages are varied on gates labeled 
c ("coupling" between dots), bp ("big dot plunger"), sp ("small dot plunger"), and n ("nose", which opens or closes the big dot). All 
transport measurements presented in Fig. |3 H] and 0are measured with source and drain connected as shown, in a magnetic field of 
130 mT normal to the plane of the 2DEG. (b) With gate voltage n = 0, the large dot opens into an infinite reservoir. Arrows mark regions 
where the small dot has an unpaired spin, leading to enhanced conductance at 12 mK (black) compared to 50 mK (red). Fitting the 
temperature dependence of the conductance (b inset), we find that Kondo temperature ranges from 110 to 300 mK: see Supplementary 
Information. In (c), the data from (b) are shown for stronger tunnel coupling to the right lead: c = — 244mV instead of — 282 mV. From 
temperature dependence of (c), we find that Tk ranges from 180 to 320 mK. 



the conductance decreases (Fig.Efb)). This, and all other measurements reported here, are performed in a magnetic 
field B = 130 mT normal to the plane of the heterostructure. The orbital effect of this modest field suppresses direct 
transmission through the small quantum dot, which we found to yield Fano lineshapes at zero magnetic field (cf. 22]). 
Due to the small g-factor of electrons in GaAs/AlGaAs heterostructures, \g\ < 0.4, the Zeeman effect of the field is 
unimportant in both 1CK and 2CK regimes - see Supplementary Information for details. 

The conductance as a function of temperature (e.g. Fig.|2b) inset) matches the expected form g(T) for a quantum 
dot in the Kondo regime, with the addition of a constant offset a: 



go < 2e 2 /h reflects the intentionally-imposed asymmetry of coupling to the two conventional leads that comprise the 
infinite reservoir (see Supplementary Information.) The normalized temperature-dependence of conductance /(T/Tk) 
is universal in the Kondo regime - it has no analytic form, but ranges from zero at high temperature to 1 at low 
temperature, with a broad logarithmic rise around T = Tk [23I 0|. A conductance offset such as we observe has 
been seen by other experimentalists j2{|, and is generically expected in the presence of potential scattering [2^. As 
expected, Tk varies strongly across the Kondo valley in both cases, and Tk is higher when the dot is more strongly 
coupled to the right lead, which increases the total T of the system. All results presented in this paper are for this 
same electron occupancy, although we have observed similar behavior in the next Kondo valley (two fewer electrons 
in small dot), as well as upon thermally cycling the device. 

Differential conductance g(T, Vd s ) — dI/dVa s is enhanced near zero bias (Fig. Etb) and (e)) when the electrostatic 
potential of the small dot is set to the middle of the Kondo valleys in Fig. H£b) or ( c )j respectively. This is a 
manifestation of the enhanced density of states at the Fermi level, widely accepted as one of the classic signatures of 
the Kondo effect, demonstrating clearly that the small dot acts as a magnetic impurity. Remarkably, the zero bias 
enhancement changes to zero bias suppression as gate n is made more negative, closing off the big dot to form a finite 
reservoir with integer occupancy (Fig. E{f)). The change signals that the single-channel Kondo state with the leads 
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FIG. 3: The formation of two competing 1CK states with two different reservoirs, (a) sp is set so that the small dot has an unpaired spin 
(middle of a Kondo valley), and c is set to — 282 mV. Conductance is shown as a function of n and bp. With n > — 180 mV, varying bp 
has little effect on the conductance. As n is made more negative, the reservoir becomes finite and a series of stripes in the conductance 
map reflects Coulomb blockade in the finite reservoir (see Supplementary Information for explanation of the modulation.) In (b) and (c), 
differential conductance as a function of bias (Va s ) an d bp is shown for n = — 170 mV and n = — 224 mV, respectively. In (b), there is 
clear zero-bias enhancement, consistent with Kondo effect. In (c), the zero-bias enhancement (denoted by righthand arrow) is modulated 
by the charge state of the finite reservoir, for reasons explained in the text [2nd paragraph after Eq. (5)]. The lefthand arrow marks 
the charge degeneracy point of the finite reservoir, (d) Same as (a), except with stronger interdot coupling: c = — 244 mV. (e) As in 
(b), we observe a clear zero-bias conductance enhancement for n = — 170 mV. (f) However, in this case at n = — 224 mV the zero-bias 
enhancement is replaced by zero-bias suppression (denoted by righthand arrow). Here the local spin forms a Kondo state with the finite 
reservoir, suppressing conductance through the small quantum dot (see Text.) Again, the zero-bias feature is modulated by the charge 
state of the finite reservoir, and the lefthand arrow marks the charge degeneracy point of the finite reservoir. 



has been broken, to form instead solely with the finite reservoir. This occurs for Jf r > J lri as shown in more detail 
in Fig. 0] below. With slightly weaker coupling to the finite reservoir (Fig. 02c), J\r > Jfr, the Kondo state is formed 
solely with the infinite reservoir. This effect requires that the finite reservoir have integer occupancy, i.e. the device 
must be set to a Coulomb blockade valley of the finite reservoir. 

In Fig. 01 we provide further evidence that, with the finite reservoir formed, two independent ICK states can exist 
depending on the relative coupling of the small dot to the two reservoirs. We have fine control over the occupancy 
of both the finite reservoir and the small dot with gates bp and sp, as shown in Fig. 0^a) and more completely in 
Supplementary Information. Conductance g{T) = g(0,T) at weak coupling to the finite reservoir (J; r > Jf r ) fits the 
expected empirical form of Eq. (3). Fig. ^Jc) shows that the conductance g(T) at many points in (sp,bp) (Fig. Ufa)) 
can be collapsed onto a universal curve. The differential conductance g(V ds , T) of a ICK system is further expected to 
follow a specific form as a function of both bias and temperatures/at an energy scale substantially below fcT K 0: 



g(0,T) - g(V ds ,T) feV dl 

— K 



T a V kT 



(4) 



where the exponent a = 2 is characteristic of ICK, and k = 0.82M-. The numerical prefactor of order unity is 
dependent on the underlying model, numerical calculations, and proximity to the symmetric 2CK fixed point (see 
Supplementary Information), so we simply treat k as a free fitting parameter for each set of gate voltages. Fig. 0Jd) 
demonstrates excellent ICK scaling at temperatures of 12, 24, 28, and 38mK, all well below Tk- A nonlinear fit to 
the data in Fig.^d) yields a = 1.72 ± 0.40 (95% CL), consistent with a = 2. 

In Fig. Efe-h), we demonstrate that at stronger coupling to the finite reservoir (Jf r > J; r ) the small dot forms a 
Kondo state with the finite reservoir, as manifested by low-energy suppression rather than enhancement of conductance 
between the normal leads of the small dot. We must modify the form we use to fit the temperature dependence to 
reflect this inversion: 

g(T) = .go (f - f(T/T K )) +a. (5) 
. ' 
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FIG. 4: Energy dependence of Kondo effect, with the finite reservoir formed, (a) - (d) Antiferromagnetic coupling to the infinite reservoir 
(normal leads) is stronger than coupling to the finite reservoir: c = -282 mV. (a) Conductance as a function of gates sp and bp, at 
Vd s =0. (b) Any point in (sp, bp) such that the occupancy of the small dot is odd and that of the finite reservoir is integer shows enhanced 
conductance at low temperature and low bias (inset: 12 mK, grey, to 60 mK), and a temperature dependence consistent with Kondo effect, 
(c) The normalized dependence of conductance on temperature is uniform for different points in (sp, bp) space, while Tk ranges from 
50mK to 180mK. g = g — a is the conductance with the temperature-independent offset subtracted off. (d) Plotting a specific combination 
of temperature and bias collapses the data for a single point in (sp,bp) space (Tk = 175 mK, go = 0.75e 2 //i) onto a single V-shaped 
curve, corresponding to the scaling relation predicted for 1CK (Eq. (e) - (h) Kondo effect with the finite reservoir: c = -244mV. (e) 

Conductance as a function of sp and bp: conductance is now suppressed rather than enhanced at low bias and temperature (cf. Fig.lSTfV) 
(f) Fitting the conductance as a function of temperature to the empirical form we expect for 1CK with the finite reservoir (Eq.[3J we find 
that Tk ranges from 30mK to 130mK. (g) We again normalize and collapse the temperature dependence at several points in (sp, bp) onto 
a single curve, (h) We collapse differential conductance data like those in inset (f) (12 mK, grey, to 30 mK) at a single point in (sp, bp) 
onto a single inverted V-shaped curve using the same temperature- bias scaling as in (d). Deviations from perfect scaling may be related 
to the slightly lower Kondo temperature (Tk = 120mK, go = 0.16e 2 //i.) 



Again temperature dependence at multiple points in (sp, bp) (Fig.^Je)) collapse onto a single normalized Kondo form 
g/go vs T/Tk (Fig. E[g)), providing strong evidence that a distinct 1CK state has formed with the finite reservoir. 
Furthermore, using the same scaling relation as above (Eq. (@J) the data again collapse onto a single (inverted) curve 
at low bias and temperature (Fig. Efh)). Interestingly, we find the numerical prefactor of n to be 0.25 in both 1CK 
regimes - precisely matching each other but only roughly agreeing with our predicted value of 0.82. 

Having established the existence of two distinct Kondo ground states - depending on the ratio Ji r /Jf r - we next 
demonstrate the tunability necessary to reach the symmetric 2CK state, J; r ~ Jp r . By setting the tunnel coupling 
to the finite reservoir to an intermediate value, we can observe either zero-bias enhancement or zero-bias suppression 
(marked by white arrows in Fig. 5(a)), in both cases away from any charge degeneracy point of the finite reservoir 
(marked by a black arrow in Fig. 5(a).) This is expected|l7l|. since the antiferromagnetic coupling to a reservoir 
depends not only on a tunneling rate but also on the energy required to transfer an electron from the local site to that 
reservoir. Gate bp tunes that addition energy for the finite reservoir, modifying Jf r while keeping J; r nearly constant. 
The region in (sp, bp) of suppressed conductance (red) grows rapidly with increasing coupling to the finite reservoir, 
as seen in Fig. OJb, c, and d) for c = —258, —256, — 254 mV, respectively. 

The evolution of dl /dVd s from zero-bias enhancement to zero-bias suppression as a function of coupling gate c may 
be seen most clearly in Fig.|3Je). We identify the curve for c = — 260 mV as being very close to the 2CK symmetric 
point. At first blush, it is surprising that this curve does not display a clear cusp at low Vd s (G ~ const — ylQ. 
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FIG. 5: Evidence for 2CK physics, (a) Differential conductance as a function of bp and Vd s in the middle of a Kondo valley, with 
intermediate coupling to the finite reservoir c = — 258 mV. In contrast to Figljjc) and|JJf), here we observe both zero-bias enhanced 
and zero-bias suppressed conductance (marked by bottom and top white arrows, respectively) by fine-tuning bp. The black arrow marks 
the charge degeneracy point of the finite reservoir, (b)-(d) At c = —258, —256, and — 254 mV, conductance may be either enhanced or 
suppressed at zero bias, depending on the fine tuning of the electrostatic potentials of the small dot and finite reservoir with gates sp 
and bp, respectively. Red and blue indicate regions of suppressed conductance or enhanced conductance, respectively, while gray indicates 
relatively flat conductance around V da = 0. White regions are too close to charge degeneracy of dot or finite reservoir to diagnose 
Kondo-induced enhancement or suppression of conductance. Increased coupling to the finite reservoir expands the region of suppressed 
conductance (red). The yellow dashed line in (b) shows the setting of sp used in (a), and the yellow dot shows the approximate location 
in the charging hexagon of the symmetric 2CK point analyzed in (f) and (g). (e) Differential conductance near zero bias evolves with 
coupling c from zero-bias enhancement to zero-bias suppression. For each curve, sp sets the small dot in the middle of a Kondo valley 
and bp sets the finite reservoir midway between two charge degeneracy points. (f,g) Tuning bp near — 304.7 mV (bottom white arrow in 
(a)), we observe that differential conductance depends on bias and temperature with a = 0.5, consistent with 2CK (f) and inconsistent 
with ICK (g, which attempts to apply scaling to exactly the same data, but with a = 2.) In the Supplementary Information we show the 
converse, namely that in the ICK region the 2CK scaling law does not fit, while the ICK scaling does. A two-dimensional nonlinear fit 
to the data set used for (f) and (g) yields «2 = 0.62 ± 0.21, consistent with a = 0.5. In (f) and (g) we do not scale the vertical axis by 
TkiT^CK, or go because here we lack independent measures of these parameters. 



In fact, conformal field theory predicts that at the 2CK symmetric point the differential conductance should depend 
quadratically on bias for eVds < kT, and should only cross over to behavior at higher bias eVds ^ 3fcT (see green 
curve in Fig- EJf ) ) ■ Such a crossover is hard to see in a single plot of differential conductance versus bias. Instead, 
we combine the dependence of differential conductance on both bias and temperature in a scaling plot to pro duce 
compelling evidence for 2CK. The expected scaling form is somewhat different from that for ICK |cl Il4l Ell Eif : 

9 (0,T)- 9 (V ds ,T) =K2Y (eV t \ (6) 



T a 2 y kT 

Here a 2 = 0.5, n 2 = {g /2){n /T 2CK ) a \ and 

\rt \ 1 n i i \ \ - 1 for X > I . . 

Y(x) = l-F 2CK (xM*{^ cx2 ^ (7) 

with c ss 0.0758 30|. As with k in the ICK analysis, in practice we treat k 2 as a free parameter for each set of gate 
voltages. F 2 ck is found by conformal field theory [8l I3lll32j . 

FigureEJf) shows that when we tune close to the 2CK symmetric point (c = —258 mV rather than —260 mV, due to 
a small shift in parameters) data at various temperatures and biases collapse onto each other and match the conformal 
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field theory prediction (Eq. (JBJ), which is scaled vertically by «2- The horizontal axis is plotted as (eV^s/kT) - 5 , to 
emphasize that g(Vd s ) ~ const — y/eVd s for eVds/kT 3> 1. 

In Fig.OJg) we show the same data scaled as would be appropriate for 1CK (Eq. |@}) instead of 2CK. As anticipated, 
this scaling fails completely: scaled data for different temperatures deviate from each other even near zero bias. A 
two-dimensional nonlinear fit to the data in Fig. [SJf) produces a value c*2 = 0.62 ± 0.21 (95% CL), consistent with 
2CK behavior. Naively, we would expect 2CK behavior to persist only up to {kT, eVds} ~ {kT^) 2 / Ec ~ 1.7/xeV [33| . 
Empirically, 2CK persists to much higher energies: conductance follows the 2CK scaling form up to Vd s = 15/ieV, 
corresponding to T = 180 mK, even higher than Tk- Enhancement of 2CK energy scales has been predicted in our 
geometry in the presence of charge fluctuationsjlij], but is not expected to be so dramatic for our parameter values. 

In this Article, we have presented data demonstrating the existence of two independent 1CK states, along with 
a study of the associated 2CK state. Remarkably, the conductance of the symmetric 2CK state matches not only 
a simple power law but rather a complete theoretically calculated non-Fermi-liquid scaling function over a broad 
range of energy (Eq. 7, Fig. 5(f)). In future, it would be interesting to extend this scaling form theoretically and 
experimentally to cover the effects of a Zeeman field and slightly asymmetric coupling to the two reservoirs, for 
example to quantitatively describe the family of curves in Figure Eff). Finally, other parameter regimes of the same 
structure may show additional exotic behavior |Ts|. 
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